Correlations in superstatistical systems 



Christian Beck 



School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London El 

4NS, UK 



Abstract. We review some of the statistical properties of higher-dimensional superstatistical 
stochastic models. As an example, we analyse the stochastic properties of a superstatistical model 
of 3-dimensional Lagrangian turbulence, and compare with experimental data. Excellent agreement 
is obtained for various measured quantities, such as acceleration probability densities, Lagrangian 
scaling exponents, correlations between acceleration components, and time decay of correlations. 
We comment on how to proceed from superstatistics to a thermodynamic description. 
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A SHORT REMINDER: WHAT IS SUPERSTATISTICS? 



Complex systems often exhibit a dynamics on two time scales: A fast one as represented 
by a given stochastic process and a slow one for the parameters of that process. As a 
very simple example consider the following linear Langevin equation 

v = -yv + oL(t) (1) 

with parameters y, o that fluctuate on a long time scale [1]. It describes the velocity v of 
a Brownian particle that moves through spatial 'cells' with different local inverse tem- 
perature /3 := y/{2o 2 ) in each cell (a nonequilibrium situation). Assume, for example, 
that the probability distribution of /3 in the various cells is a ^-distribution of degree n, 

/(/3)~/3' 3 / 2 -V*. (2) 

Then the conditional probability is p(v\f$)) ~ e~^ v2 , the joint probability is p(v, j8) = 
/(/3)p(v|j8), and the marginal probability is p(v) = J ° o /(j8)p(v|/3) J/3. Integration yields 

p(v) ~ o+yito-iy)'/^) <3) 

i.e. we obtain power-law stationary distributions just as in Tsallis statistics [2] with 
q =l + -L T J= 2p /(3-q), where #> = //(jB)j3</jB is the average of j8. 

All this has a very broad interpretion and can be generalized in various ways — /3 
need not to be inverse temperature. One can generalize the above example to general 
probability densities /(j8) and general Hamiltonians in statistical mechanics. One then 
has a superposition of two different statistics: that of /3 and that of ordinary statistical 



mechanics. The short name for this is superstatistics yfl. Superstatistics describes com- 
plex nonequilibrium systems with spatio-temporal fluctuations of an intensive parameter 
(e.g. inverse temperature) on a large scale. 
Define an effective Boltzmann factor B(E) by 

B(E) = / f(P)e^ E dP 



where /(j8) is the probability distribution of /3 and E the energy of the system. Many re- 
sults can be proved for general /(/3) . Here we list some recent theoretical developments 
of the superstatistics concept: 

• Can prove superstatistical generalizations of fluctuation theorems [0] 

• Can develop a variational principle for the large-energy asymptotics of general 
superstatistics [5] (depending on /(/3), one can get not only power laws for large E 
but e.g. also stretched exponentials) 

• Can formally define generalized entropies for general superstatistics [@, @] 
Can s tudy va rious theoretical extensions and workouts of the superstatistics concept 

HI EES 



• Can prove a superstatistical version of a Central Limit Theorem leading to Tsallis 
statistics 

• Can relate it to fractional reaction equations [fl5ll 

• Can consider superstatistical random matrix theory lfl6ll 



• Can apply superstatistical techniques to networks [|17|1 and time series [|18|1 
...and some more practical applications: 

• Can apply superstatistical methods to analyse the statistics of 3d hydrodynamic 
turbulence flM [H H IzE S 



Can apply it to atmospheric turbulence (wind velocity fluctuations at Florence 
airport JlElll]) and defect turbulence Q 
Can apply superstatistical methods to finance | 26l l27| 
Can apply it to solar flares ll28ll . and even to print queues [H^ 
Can apply it to cosmic ray statistics [(3( 



Can apply it to various scattering processes in particle physics ||3lL l32| l 
Can apply it to hydroclimatic fluctuations ||3~3"|l 



Can apply it to British train delay statistics D34Q 



PHYSICALLY RELEVANT SUPERSTATISTICAL 
UNIVERSALITY CLASSES 



Basically, there are 3 physically relevant universality classes [18]: 

• ( a ) Z 2 -superstatistics (= Tsallis statistics) 

• (b) inverse ^-superstatistics 



• (c) lognormal superstatistics 



Why? Consider, e.g., case (a). Assume there are many microscopic random variables 
j = 1, ... ,7, contributing to /3 in an additive way. For large 7, their sum -jjY?j=\ %j 
will approach a Gaussian random variable X\ due to the (ordinary) Central Limit The- 
orem. There can be n Gaussian random variables X\, . . . ,X n due to various relevant de- 
grees of freedom in the complex system. Since /3 is positive we may square the X,- to 
obtain something positive. The sum /3 = YJi=\Xf ls then ^-distributed with degree n, 
i.e., 

1 ( n \ n l 2 , "P 

/G6) = f4t j3 M/2 -V4, (4) 



r(|) V2/3o 

where /3o is the average of /3. Integration as described in section 1 yields Tsallis statistics 
as a special case of superstatistics. 

(b) The same considerations can be applied if the 'temperature' j6 rather than /3 
itself is the sum of several squared Gaussian random variables arising out of many 
microscopic degrees of freedom The resulting /(j8) is the inverse ^-distribution: 

go ( "fc\" /2 r „/2-2 



/(j8) = f% \ 2 ) r 6 ~ W - (5) 

It generates superstatistical distributions p(E) ~ / f(fi)e~P E that decay as e~$"^ for 
large £ |5|]. 

(c) /3 may be generated by multiplicative random processes. Consider a local cascade 
random variable X\ = IT/=i <■>./> where 7 is the number of cascade steps and the E,j are 
positive microscopic random variables. By the Central Limit Theorem, -^logXi = 
^£j =1 log^j becomes Gaussian for large 7. Hence X\ is log-normally distributed. 

In general there may be n such product contributions to /3, i.e., /3 = rT/Li^'- Then 
log/3 = Y!i=\ l°gX; is a sum of Gaussian random variables; hence it is Gaussian as well. 
Thus /3 is log-normally distributed, i.e., 

/(/3) = ^=L^exp<|^f^ j>. (6) 



Lognormal superstatistics is relevant in turbulence |18, T^, 20. 21. 221 



APPLICATION TO LAGRANGIAN TURBULENCE 



Turbulence is a spatio-temporal chaotic state of the Navier-Stokes equation. Energy is 
dissipated in a cascade-like process. Bodenschatz et al. ||22l. B5L |36|] obtained rather 
precise measurements of the acceleration a(t) of a single tracer particle in a turbulent 
flow. One can now construct a superstatistical Lagrangian model for 3 -dimensional 
velocity differences u(t) :=v(t + t) — v(t) of such a tracer particle (note that a = u/x for 



LN-5uperata1is1ics' 

'R-200' + 

■R-CSO 1 X 

'R-970' * 




20 40 



FIGURE 1. Distribution of acceleration as measured by Bodenschatz [2 
eq. ®, s 2 fa 3 . 



and as predicted by 



small t). This model is given by the superstatistical stochastic differential equation 1190 



-yu + Bn x u + oL(t). 



(7) 



The new thing as compared to previous work is the term involving the vector product. 
It describes fluctuating enstrophy (rotational energy) around the test particle. While y 
and B are constants, the noise strength o and the unit vector n evolve stochastically on a 
large time scale T a and 7g, respectively. One has T a y^R^ >> 1 ? where is the Taylor 
scale Reynolds number. The time scale describes the average life time of a region of 
given vorticity surrounding the test particle. 

Define /3 := 2y/a 2 , then in this model /3 1 ~ v 1 / 2 (e) _1//2 e, where v is the kinematic 
viscosity and (e) the average energy dissipation. The probability density of the stochastic 
process /3 (t) is assumed to be a lognormal distribution as given in eq. ©. For very small 
T an acceleration component of the particle is given by a x = u x J % and one gets the 
following prediction for the stationary distribution: 



P(a x ) 



2ns 



dp j8 



-1/2 



exp 



2s 2 



P\2 



(8) 



This compares very well with the experimentally measured probability distribution of 
acceleration, see Fig. 1. 



CORRELATIONS INDUCED BY SUPERSTATISTICS 



3-dimensional superstatistics induces correlations between the components a x ,a y ,a z of 
the acceleration vector a. Consider the ratio R : = p(a Xl a y ) / (p(a x )p(a y )). For indepen- 
dent acceleration components this ratio would always be given by R = 1 . However, our 



FIGURE 2. The quantitity R as given by eq. (O for lognormal superstatistics. 



3-d superstatistical model yields the prediction 

This is a very general formula, it is also valid for Tsallis statistics, where /(/3) is the 
^-distribution. Note that R = 1 for /(j8) = <5(/3 — j8o), i.e. if there are no fluctuations in 
/3 then all components are independent random variables. 

Fig. 2 shows R = p(a x ,a y ) / (p(a x )p(a y )) as predicted by lognormal superstatistics. 



The figure strongly resembles R as experimentally measured by Bodenschatz et al. 112211 
in a turbulent flow, see Fig. 3. 

Besides correlations between components one can also look at temporal correlations. 
The superstatistical model tl9fl allows for the calculation of temporal correlation func- 
tions as well. In particular, we may be interested in temporal correlation functions of 
single components u x of velocity differences, i.e. C(t) = (u x (t r + t)u x (t')). By averaging 
over the possible random vectors n one arrives at the formula 

C(t) = ^(^)e _1 *(2cosft+l), (10) 

i.e. there is rapid (exponential) decay with a zero-crossing at t* = Exponential 
decay and zero-crossings are also observed for the experimental data. The model 
also correctly reproduces the experimentally observed fact that the correlation function 
of the absolute value \a\ decays very slowly as compared to that of the single compo- 
nents. Moreover, it correctly describes the fact that enstrophy lags behind dissipation 
& 




FIGURE 3. iog lQ R as measured by Bodenschatz et al. Iq^R . 



LAGRANGIAN SCALING EXPONENTS 

The moments of velocity differences of a single Lagrangian test particle that is embed- 
ded in a turbulent flow scale differently from those measured in a fixed laboratory frame. 
Our superstatistical model lfl9ll allows for the analytic evaluation of the Lagrangian scal- 
ing exponents. The moments of velocity difference components u x on a time scale T are 
obtained as 

( U J X ) = (j-l)\\ m -iw* j2 . (11) 

Assuming simple scaling laws of the form m ~ x a , w = e sl ~ x b , where a and b are 
so far arbitrary real numbers, on gets (u{) ~ ~ T~ a ? +b %J a . Hence the Lagrangian 
scaling exponents are given by £ 7 - = — |j + |j 2 . Usually one assumes ^2 = L hence we 
get a = jb — 1 thus 

Zj = {\+X 2 )j-\x 2 j\ (12) 

where A 2 := — \b. This prediction is in good agreement with the recent measurements 
of Bodenschatz et al. [36], see Fig. 4. 

FROM SUPERSTATISTICS TO (GENERALIZED) 
THERMODYNAMICS 

We end this paper with some more general thoughts. Can we proceed from superstatistics 
as a merely statistical technique to a proper thermodynamic description! There are some 
early attempts in this direction by Tsallis and Souza §@\. Here we want to follow a 
somewhat different approach [7]: One starts quite generally from two random variables 
E and B (representing energy and invere temperature) and then considers the following 
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FIGURE 4. Lagrangian scaling exponents. Data points: Measurements of Bodenschatz et al. [ 36] . Solid 
line: Theoretical prediction of the superstatistical model (A 2 — 0.085) lfl9ll . Dashed lines: Some other 
competing models 113811 . 



effective entropy for a superstatistical system 
S[E,B] = S[E\B]+S[B] 

= J dPf(P)(pu(P)+\nZ(P))-jdPf(P)\nf(P) 1 

where U is the local internal energy and Z the local partition function. One can do ther- 
modynamics with this extended entropy function. It reduces to ordinary thermodynamics 
for f(P) = 8(j3 — fio). For sharply peaked distributions /(/3) this is a slightly deformed 
thermodynamics, which can be evaluated in a perturbative way. One can also maximize 
this entropy with respect to appropriate constraints in /3 to get e.g. a lognormal distribu- 
tion for f(P), or generally some other distribution /(/3) depending on the constraints. 
For more details, see OD. 

SUMMARY 

• Superstatistics (a 'statistics of a statistics') provides a physical reason why more 
general types of Boltzmann factors (e.g. of power-law type) are relevant for 
nonequilibrium systems with fluctuations of an intensive parameter. 

• There is evidence for three major physically relevant universality classes: % 2 - 
superstatistics = Tsallis statistics, inverse # 2 -superstatistics, and lognormal super- 
statistics. These arise as universal limit statistics for many different complex sys- 
tems. 

• Superstatistical techniques have been successfully applied to a variety of complex 
systems. 




• A superstatistical model of Lagrangian turbulence |fl9ll is in excellent agreement 
with the experimental data for probability densities, correlations between compo- 
nents, decay of correlations, and Lagrangian scaling exponents. 

• The long-term aim is to find a good thermodynamic description for general super- 
statistical systems. 
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